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We study the generation of primordial perturbations in a (single-field) slow-roll 
inflationary universe. In momentum space, these (Gaussian) perturbations are char- 
acterized by a zero mean and a non-zero variance A'^{k,t). However, in position 
space the variance diverges in the ultraviolet. The requirement of a finite variance in 
position space forces one to regularize A^(fc, i). This can (and should) be achieved 
by proper renormalization in an expanding universe in a unique way. This affects 
the predicted scalar and tensorial power spectra (evaluated when the modes ac- 
quire classical properties) for wavelengths that today are at observable scales. As 
a consequence, the imprint of slow-roll inflation on the CMB anisotropics is sig- 
nificantly altered. We find a non-trivial change in the consistency condition that 
relates the tensor-to-scalar ratio r to the spectral indices. For instance, an exact 
scale-invariant tensorial power spectrum, rit = 0, is now compatible with a non-zero 
ratio r ~ 0.12 it 0.06, which is forbidden by the standard prediction (r = —8nt). The 
influence of relic gravitational waves on the CMB may soon come within the range 
of planned measurements, offering a non-trivial test of the new predictions. 

PACS numbers: 98.80.Cq, 04.62. +v, 98.70.Vc 



I. INTRODUCTION AND SUMMARY 

Inflation [1] provides a natural solution to the horizon and flatness problems of the hot 
big-bang cosmology. A sufficiently long period of rapid expansion can explain the large 
scale homogeneity, isotropy, and flatness of our visible universe. Inflation also provides 
a quantitative explanation [2J to account for the origin of small inhomogeneities in the 
early universe. These inhomogeneities are responsible for the structure formation in the 
universe and for the anisotropies present in the cosmic microwave background (CMB), which 
were flrst detected by the COBE satellite and further analyzed by the WMAP satellite 
The potential-energy density of the inflaton fleld is assumed to cause the inflationary 
accelerated expansion, and the ampliflcation of its quantum fluctuations and those of the 
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metric are inevitable consequences in an expanding universe j^. These fluctuations acquire 
classical properties in the inflationary period and provide the initial conditions for classical 
cosmological perturbations after the big-bang. The detection of the effects of primordial 
tensorial metric fluctuations (gravitational waves) in future high-precision measurements of 
the CMB anisotropics will serve as a highly non-trivial test of the inflationary paradigm 
and to constrain speciflc models. Therefore, it is particularly important to scrutinize the 
predictions of inflation for the tensorial and scalar power spectra. In this respect, it was 
pointed out in [5,] (see also [61]) that quantum fleld renormalization signiflcantly modifles the 
amplitude of quantum fluctuations, and hence the corresponding power spectra, in de Sitter 
inflation. The analysis was further improved in 0] (see also the essay [sl) to understand how 
the basic testable predictions of (single-fleld) slow-roll inflation could be affected by quantum 
fleld renormalization. In this work we further study this issue, improve the technical analysis, 
and provide a more complete and robust discussion of how the observable consequences of 
inflation are altered when quantum fleld renormalization is taken into account. 

Let us briefly summarize the logic of our approach. Let us assume that (p{x, t) represents 
a perturbation obeying a free fleld wave-equation on the inflationary background ds"^ = 
—dt^ + a?{t)dx^ , where a{t) is a quasi-exponential expansion factor (a(t) ~ e^^ ). At the 
quantum level, this fleld is expanded as 

¥.(f , t) = I d'k[Mt)a^k^''' + ^male-^''] , (1.1) 

where the creation and annihilation operators satisfy the canonical commutation relation 
[ag, at] = 6^{k — k'). The mode functions (fk{t) are required to satisfy the adiabatic condition 
(see, for instance, |9|). The power spectrum for this perturbation, Ai^ (fc, t), is usually defined 
in terms of the Fourier transform of the variance of the field |lOl . Ill| 

(vdmlXt)) = S'{k - k')—Al{k,t) , (1.2) 

where 0j^{t) = ipk{t)ai^. These modes describe a perturbation field characterized, in momen- 
tum space, by a zero mean {ff^it)) = and the variance (11.21) . The advantage of working 
in momentum space resides in the fact that different modes fluctuate independently of each 
other, as explicitly displayed by the presence of the delta function in (II. 2p . This way, the 
quantum fleld is regarded as an inflnite collection of oscillators, each with a different value 
of k. In position space the perturbation is also characterized by a zero mean {(p{x,t)) = 
and a variance (or dispersion) 

{^'{x,t)) = j^,J rf^fcd=^fc'(^,Kt)4,(i))e^^'^-'''^" , (1-3) 

which, due to spatial homogeneity, turns out to be independent of x. This variance is 
formally related to the power spectrum by 

dh 

yAj(A;,t). (1.4) 

As is well-known in quantum fleld theory, the above expectation value quadratic in the fleld 
(p is divergent. It suffers from quadratic and logarithmic ultraviolet divergences 

{V>\x,t)) — dk{- + \ + ...). (1.5) 

47r^ Jo a^k 
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The first term corresponds to the usual contribution from vacuum fluctuations in Minkowski 
space. This contribution can be eliminated by renormalization as claimed, for instance, in 
T^l- However, the second term is also characteristic of vacuum fluctuations in a curved 



background. Because the different fc-modes fluctuate independently of each other, one could 
be tempted to get rid of this logarithmic ultraviolet divergence by simply eliminating the 
modes with k > aH and leaving the rest unaffected (see, for instance, l3 ). If one eliminates 



this divergence using a window function in this way, as is usual for random fields, then one 
obtains A^(A;) ^ H'^/in'^, where A^(/c) is defined by the quantity A'^{k,t) evaluated a few 
Hubble times after the "horizon crossing time" t^, since this is the time scale at which the 
modes behave as classical perturbations.^ However, one should take into account that the 
field fluctuations are quantum in nature and, therefore, one should consider the subtle points 
of quantum field theory (QFT) regarding the ultraviolet divergences. 

Even though free quantum field theory is usually regarded as an infinite set of independent 
harmonic oscillators (one for each k-mode), there are fundamental holistic aspects of QFT 
that can not be properly understood in terms of independent modes. Renormalization is 
the hallmark of the holistic aspect of QFT. This is clear in the fact that, although the 
renormalization schemes in QFT in curved spacetimes are based on the ultraviolet behavior 
of the theory, the infrared sector is also affected by renormalization, leading potentially to 
observable consequences. This can be explicitly displayed by considering, for instance, the 
Casimir effect. The energy density between the two conducting plates obtained by proper 
renormalization provides the well-known and experimentally tested expression. However, 
a naive subtraction obtained by introducing a high-frequency cut-off in the integrals in 
momentum space (i.e., treating the fc-modes as being independent) produces a quite different 
result (see the discussion of jl4|). 

Taking this into account, we see that the logarithmic divergence in (11.51) should be dealt 
with by renormalization and one can not rule out, a priori, the possibility that the treatment 
of the divergences at very high values of k may produce some impact at lower momenta. 
Therefore, we propose that in the standard definitions of the spectrum A^(fc,t), as given 
in (11.21) and (II. 4p . one should replace the unrenormalized {(p'^{x,t)) by the renormalized 
variance, {ip'^{x,t))ren- Writing A'^{k,t) for the spectrum defined in this way, the definition 
in (11.41) (and similarly in (11.21) ) is replaced by the corresponding renormalized expression 







dk ~ 

{if\x,t))ren= / yAj(fc,t) . (1.6) 



This way, the physical variance {(p'^{x,t))ren remains a well-defined quantity, in the same 
way as one could obtain a finite expression for the expectation values of the quantum stress- 
energy tensor. To complete the physical consistency of this approach, it would be desirable 
to define a unique expression for the necessary subtractions required to produce a consistent 
{ip'^{x,t))ren- Slucc the power spectrum is defined in momentum space, the natural scheme 
is renormalization in momentum space, so we define 

4-Tr P'^ 



{^\x,t))ren = J^J^ k' dk{\MW ' C,it)) , (1.7) 



^ The time ifc is defined by a{tk)/k = H{tk), where a{t) is the expansion factor and H ~ d/a is the Hubble 
rate. 
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where OJt) represents the expected counterterms. As we will see, adiabatic renormalization 
[sl, E, [3] , which works by subtracting a set of counterterms "mode- by-mode" , provides a 
natural expression for the counterterms encoded in Ck{t). Moreover, the DeWitt-Schwinger 
renormalization, originally defined in position space (see, for instance, [ol, 16]), can be nicely 
translated to momentum space thus providing another answer for Ckit). When these 
two schemes are applied to the field perturbations arising from inflation, the resulting ex- 
pressions for Ck{t) coincide, thus defining a unique expression for the spectrum A^(A;,t). 

The holistic nature of QFT is then explicitly realized through (11 .Op . Although the coun- 
terterms are fully determined by the ultraviolet behavior of the modes, the long wavelength 
sector, and hence the new A^(A;, t), is significantly affected by the subtractions. In the 
slow-roll scenario, when H slowly decreases with time, the effects of renormalization have a 
non-trivial impact on A^(fc,t) when this quantity is evaluated a few Hubble times after the 
time of horizon crossing tk- For instance, for the tensorial modes we obtain 

where e is the usual slow-roll parameter. A similar expression is obtained for the scalar 
perturbations (involving now the slow- roll e and rj parameters), with the corresponding 
changes in the tensor-to-scalar ratio, the spectral indices and the consistency relation. The 
new predictions remain in agreement with observation for the simplest forms of inflation 
(0^ and 0^ potentials). It is worth pointing out that the new consistency condition for 
single-field inflation is predicted to be 



An' 



4(1 -n,- nt) + ^ [l-Us- sj2n[ + (1 - n^f - ngj , (1.9) 

instead of the standard prediction r = — Sn^. The tensor-to-scalar ratio is now related 
with the spectral indices nj, (1 — ng), and also n[ = drit/dlnk. This modification has far 
reaching consequences. For instance, since the observations from the 5- year WMAP [3'](with 
BAO+SN) strongly suggest that (1 - n,) ^ 0.030±0.015 (with r < 0.22), expression (K^ 
allows for an exact scale invariant tensorial power spectrum, nt = 0, while being compatible 
with a non-zero ratio r ^ 0.12 ± 0.06. We will comment on this further later on. 

The paper is organized as follows. In section II we briefly review the standard ways 
to derive the expressions for the tensor and scalar power spectra of single-field slow-roll 
infiation. In section III we work out the new definition of the power spectra and give the 
technical details involving the required renormalization in momentum space. In section 
IV we provide, as a consequence of the new power spectra, the corresponding expression 
for the tensor-to-scalar ratio r, the tensorial and scalar spectral indices, and the slow- roll 
parameters. This leads, in particular, to a change in the consistency condition that relates 
the tensor-to-scalar amplitude ratio to the spectral indices. In section V we summarize our 
results and conclusions. We use natural units h = 1 = c. 



II. SPECTRUM OF FLUCTUATIONS FROM INFLATION 

We will now proceed to briefiy review the standard predictions for the power spectra in 
single-field slow-roll infiation. In obtaining these predictions, quantum field renormalization 
in the curved spacetime of the expanding universe is not taken into account. 
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A. Tensorial spectrum 

Let us focus on the production of rehc gravitational waves by considering fluctuating 
tensorial modes hij [x, t) in an expanding, spatially fiat universe 

ds^ = -dt^ + a^{t){5ij + hij)dx'dx^ . (2.1) 

The wave equation obeyed by these modes comes from the linearized Einstein equations and 
is given by 

— a^hij — Sadhij + V^hij = . (2-2) 

Expanding the fluctuating fields h^j in plane wave modes hk{t)eije^^^ , where Cjj is a constant 
polarization tensor obeying the conditions Cjj = Cji^eu = and kiCij = 0, we obtain the 
equation 

hk + 3Hhk + = , (2.3) 

with k = \k\. The conditions for the polarization tensor imply that the perturbation field 
hij can be decomposed into two polarization states described by a couple of massless scalar 
fields h^^y.{x,t), both obeying the wave equation fl2.3p [l8|] (see also, for instance, 19 



from now on, we omit the subindex + or x). On scales larger than the Hubble radius, 
the damping term 3Hhk dominates. However, on scales smaller than the Hubble radius 
it is the spatial gradient term that dominates over the damping term, thus leading to the 
conventional flat-space oscillatory behavior of modes. To constrain the form of the modes 
defining the quantization, it is natural to impose the adiabatic asymptotic condition 0, 0] 
for large k 

~ (2(27r)=^w(t)a3(t))-V2e-^r-(*')'i*' ^ (2.4) 

V 167rG 

with w{t) = k/a{t), where the factor vTGttG {G is the Newton constant) has to be introduced 
to get a canonically normalized variable. This condition does not uniquely fix the form of 
the modes. For instance, in an exact de Sitter background {a{t) = c^*, with H a stric 
constant) one can invoke de Sitter invariance to uniquely determine the modes [20|, and one 
obtains 

These modes oscillate until the physical wave length reaches the Hubble length H~^. The 
amplitude of the modes at this time, usually called the "horizon exit" time tk, defined by 
k/ci{tk) = H, is then = ^^r- A few Hubble times after horizon exit, the modes 



get frozen as classical perturbations |10l . Il3l . Il9l | (see also |2l|) with constant amplitude 



l^fcP = The freezing amplitude is usually codified through the unrenormalized quantity 
A\{k) = 47r/c^|/ifcp. Taking into account the two polarizations, one easily gets the standard 
scale free tensorial power spectrum Pt{k) = 4A^(/c) = where Mp = l/y/SnG 

is the reduced Planck mass in natural units. The appearance of such frozen fluctuations 
converts the modes into classical perturbations with wavelengths that are still stretched by 
the rapid expansion to reach astronomical scales. This is essentially what happens in the 
inflationary era. 
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1. Slow-roll inflation 

To take into account that inflation lasts for a finite period of time, one usually considers 
the so-called slow roll scenario Il9| . The homogeneous part of the inflaton field 0o (^) 



rolls slowly down its potential V{(j)) towards a minimum. Both 0o and H y ^y^V^(0o) are 
changing very gradually and this change is parameterized by the slow-roll parameters e, rj, 
where e = —H/H"^, and rj — e = (f)o/ {Hcf)^). In the slow-roll approximation, defined when 
the parameters are small e, jr^l <^ 1, one can relate them to the derivatives of the inflaton 
potential e = (M|,/2)(y/V)^, = M'^iV" /V). To generate an approximate form for the 
modes, it is convenient to introduce the conformal time variable t = J dt/a{t). In terms of 
r the wave equation turns out to be of the form 

(Phu dhu , r, , ^ . 

—4 + '2Ha-^ + = , 2.6 

and taking into account that in the slow-roll approximation 

(1-.)--^. (2.7) 

we get 

^-'-^P + eh, = 0. (2.8) 

dr^ T dr 

Within this approximation, and treating now the parameter e as a constant, one can exactly 
solve the above equation as follows 

hk{t) = {-167rGTTT/A{27rfa^y/^[E{k)Hl^\~kT) + F{k)Hl^\-kT)] , (2.9) 

where the index of the Bessel function is = 3/2 + e, and the complex coefficients E{k) and 
F{k) obey the normalization requirement 

\E{k)\'-\F{k)\' = l. (2.10) 

The adiabaticity condition (12.41) implies that 

lim E{k) = 1 lim F{k) = . (2.11) 

fc— >oo fe— >oo 

The simplest way to choose E{k) and F{k) would be to require that, for e — 0, we recover 
the exact de Sitter form of the modes. This would mean that E{k) = 1 and F{k) = for 
every value of k. 



2. Infrared divergences 

The choice E{k) = 1 and F{k) = is, however, an idealized situation which assumes that 
inflation started at an infinite time in the past. A consequence of this assumption is that 
the two-point function 

{0\h{x)h{x')\0) = [ d^ke'^^'^-'''^hk{T)hl{T') (2.12) 
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is ill-defined due to an infrared divergence. This is so because in the limit — > the 
integrand in fl2.12p behaves as 

dkl 

and e > 0. For general E{k) and F{k) we have, instead 

dk\E{k)- F{k)\'^ 



(2-13) 



k 



2e 



(2.14) 



The infrared divergence is avoided if \E{k) — F{k)\'^ ^ 0, as /c — 0. This happens naturally if 
one assumes that inflation started smoothly at some early, but finite, time (see, for instance, 
(2^). Assuming that the initial vacuum is well-defined, and since the dynamical evolution 
cannot generate infrared divergences [2^, 2^, one should get E{k) F{k), when — >■ 



(as explicitly obtained in j22|), and a finite contribution to the two-point function in the 
infrared end. 

The scale of this infrared behavior of the functions E{k) and F{k) is given by the Hubble 
radius at the beginning of inflation. Only those wavelengths that were already outside the 
inflationary Hubble radius 

k/a{t) < H{t) (2.15) 

when inflation began can see this infrared behavior. However, only the wavelengths that 
crossed the Hubble radius 

k/a{tk)=H{tk) (2.16) 

at 50 — 60 e-folds before the end of the inflationary epoch have cosmological observable 
consequences today. For these wavelengths the functions E{k) and F{k) behave as in the 
Bunch-Davies vacuum {E{k) ~ l,F{k) ~ 0). Despite having no direct observable conse- 
quences, the behavior of E{k) and F{k) for very large wavelengths is important for the 
consistency of the theoretical framework because it avoids an infrared divergence in the 
two-point function. 



3. Nearly scale-invariant spectrum 

The Hubble exit time tk can be rewritten, in terms of the conformal time, as 

-kT{tk) = (1 + e) . (2.17) 

Therefore, at the Hubble exit time tk the amplitude of the modes hk given in fl2.9p takes 
the value = ^^^^x^? and then A|(fc) = §^iILJhl_ This provides a nearly scale- invariant 
spectrum of fluctuations. The k dependence of H'^{tk) 

(2.18) 



dink 

leads to 

Alik) = Aliko) i^-^J , (2.19) 

where fco is a pivot scale. One defines the tensorial spectral index rit as the exponent in the 
above expression. So 

rit = -2e{tk) . (2.20) 
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In the latter formula we have explicitly considered that the constant parameter e in the 
previous calculation is given by its value around the Hubble exit for the k mode. 

The same conclusion can be achieved by evaluating the power spectrum a few e-foldings 
after crossing the Hubble radius. One can then consider the approximation —kr -C 1, which 
implies 

\Hi'\-kT)\' -i-kr)-"' . (2.21) 

71 

Furthermore, from the slow-roll equation (12. 7p we obtain 

a = a(ro) (^^^ , (2.22) 
where tq is an arbitrary reference time. Therefore, we can write 

which explicitly shows the dependence on fc"^*^ of A^(A;) oc k^\hk{t)\'^. Now, taking for 
convenience the reference time tq in (12.221) as tq = T{tk) and using that a{tk) = k/H(tk) and 
T{tk) = — ; one can find immediately that 



\Mt)\' - :;^.H\tu) [j-^ j ^ :;^.H\t,) , (2.24) 

Taking into account the two polarizations one finally gets the same unrenormalized power 
spectrum Pt{k) = 4A^(A;) = 167rA;^|/ifcp (up to a factor 1/2) 



B. Scalar spectrum 



We now proceed to reproduce the standard results for the unrenormalized power spectrum 
of scalar perturbations. We first consider, for simplicity, the inflaton field in an unperturbed 
background and compute the evaluated at tk- In order to evaluate A^ at later times (a 
few e-folds after t^), we must improve the calculation by taking into account the fluctuation 
of the background metric using the spatially-fiat slicing. We also show how to reproduce 
this result using a more rigorous approach in terms of gauge invariant quantities. 

Consider the inflaton field (f){t, x) made out of a homogeneous part 0o(i) and a small fluc- 
tuating part 6(f){x,t). The fluctuation 6(j){x,t) satisfies, in the approximation of considering 
an unperturbed background spacetime, the wave equation 

(50 + 3H6(f) - a-^V^6(j) + V"{(j)o)6(f) = . (2.26) 

Due to the slow-rolling of (f)o{t) from the hill potential, the term V" is very small, which allows 
one to estimate the amplitude of quantum fluctuations A| in a way similar to A^. Inflaton 
fluctuations translate into curvature perturbations, which constitute the "seeds" for struc- 
ture formation and are characterized by their scalar power spectrum Pnik:) = (H/cpQ ) A?. 
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Evaluating and yH/(poj at the time of horizon exit 13, 1^ one gets = H'^(tj^)/2Tc'^ 

and H\tk)/4>o\tk) = ^M^y Therefore, Pn{k) ~ (^^) . 

One can improve this approximation by taking into account the effect of the metric per- 
turbation in ( ]2.26p . To properly define the wave equation for 6(f){x,t) we have to specify a 
slicing. Choosing for convenience the spatially-flat slicing, the equation for the field pertur- 
bation 6(j){x,t) is (see, for instance, 13|) 



S<j) + 3HS(f) - a-^V^6(j) + V"{<j)o)S(f) + \^ ( ] 6<f) = . (2.27) 

a-^ at \ H I 

In the slow-roll approximation, the corresponding Fourier components of S(f) obey 

+ 3H5'(t)k + a~^k'^S(j)k + H^i^V - 6e)50fc = . (2.28) 

To deal with this equation, one usually considers that the e and t] parameters are constant, 
while H = a/a remains a time-dependent function. So we have 

• 2 

+ 3-50fc + a-^P6(f)k + ^{Sr] - 6e)6(f)k = . (2.29) 
a 

With this assumption one can exactly solve the equation in terms of the conformal time 

r. The form of the modes 6(j)k is the same as for the tensorial ones, up to the coefficient 
VlGnG, 

6Mt) = {-rn /A(2ti fay^Hl^\-kr) , (2.30) 

where the index of the Bessel function is now = 3/2-|-3e — Evaluating this perturbation 
a few Hubble times after the horizon crossing time, one obtains the standard result 

PAk) ^ {Hi^.)\.em\^ ^ (^)^ . (2.31) 

where the constant parameter e is given at the Hubble radius time scale tk- Due to the 
fc-dependence of H{tk) and e(tfc), one also gets a nearly scale-free spectrum 

k ^ 



Pn{k) = Pniko) (2.32) 

where the so-called scalar spectral index Ug is found to be 

ns = l-6e + 2r] . (2.33) 

To complete our presentation, we will explain how the above result can be recovered 
without specifying any particular slicing and threading of spacetime. These scalar perturba- 
tions are commonly studied through the gauge-invariant quantity TZ (the comoving curvature 
perturbation) 

7^ = ^ + -^S(f) , (2.34) 

00 
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where \I/ is the curvature perturbation (i?*^^^ = 4V^\l//a^) of the spatial metric gij = — 
2'^)6ij + 2dijE]. In momentum space, TZ obeys the equation 25|, l26[ 

where z = acpo/H. In the slow roll approximation z~^dz/dT = aH{l + 2e — 7]) and hence 

-l-3e+r) 

Jo 



4r) = z{to) {-] , (2.36) 



where tq is again a reference instant of time. Furthermore, the above wave equation simplifies 
to 

d^TZk 2(1 + 3e - r?) rf7^fc , 2_ , , 

-nr-- '-^—^ + k^nk = 0. 2.37 

dT'^ T dr 

The solution obeying the adiabatic condition (and the de Sitter symmetry for H constant) 
is 

Uuit) = {-TTT/A{2T,fz^f'^H^^\-Tk) , (2.38) 

where /i = 3/2 + 3e— Evaluating TZk{t) at t = one obtains the power spectrum estimated 
above in terms of the fluctuations of the infiaton field. We can also evaluate the amplitude 
|7^fcp a few Hubble times after tk- Using r(tfc) as reference time in ( 12.361) . it is easy to get 
{Pn{k) = Al{k)=47,k'\nk\') 



C. Tensor-to-scalar amplitude ratio and spectral indices 

We have seen that the concrete value of the amplitude of the power spectra depends on 
the time at which they are evaluated. In the literature, one find different times that can be 
parameterized by A, according to the condition k/ajt) = XH. The power spectra are then 
modified by the factor (1 + A^) (see for instance |27|), which goes to unit exponentially fast. 
This ambiguity, however, is irrelevant in the evaluation of physical observables such as the 
tensor-to-scalar ratio r and the spectral indices. In fact, irrespective of the evaluation time 
{tk or a few e- foldings after it), the ratio 

Ptik) 

gives a constant quantity 

r = 16e(4) . (2.41) 

The other two physical observables are neither affected by the ambiguity in the evaluation 
time of the power spectra. The spectral indices are unambiguously given by 

nt = -2e{tk) , (2.42) 



and 



l-n, = 6e{tk) - 2r]{tk) 



(2.43) 
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From the above formulas one infers a necessary relation between measurable quantities 

r = -Srit , (2.44) 

which, as claimed many times in the literature, should be verified by any single-field slow- 
roll inflationary model irrespective of the form of the potential. For this reason, the future 
experimental checking of this condition is usually regarded as an important test of the 
simplest forms of inflation. 

This concludes our review of the standard derivation of the unrenormalized predictions 
of single-field slow-roll inflation. 



III. THE ROLE OF ULTRAVIOLET DIVERGENCES: RENORMALIZATION 
OF THE SPECTRUM OF FLUCTUATIONS 

In this section we reexamine the derivation of the power spectra taking into account the 
effects of renormalization. 



A. Ultraviolet divergences and momentum-space renormalization 



It is easy to see that Al{k,t) gives the formal contribution, per dink, to the variance of 
the gravitational wave fields h+ x 



k'dk 



dVL\hi.\'^ 



dk 
T 



Al{k,t) 



(3.1) 



since, as already defined, A\{k,t) = 4:nk^\hk\'^. Although, as noted previously, the small k 
behavior of the functions E{k) and F{k) cures the potential infrared divergence in the above 
integral, the large k behavior of the modes makes the integral divergent 



dk IGTiGk^ 
k 



I (2 + 36) 
k^ 2kH^ 



(3.2) 



It is a common view to bypass this point by regarding h{x,t) as a classical random field. 
One then introduces a window function W{kR) in the integral to smooth out the field on a 
certain scale R and to remove the Fourier modes with k~^ < R. However, as explained in the 
introduction, it is our view to regard the variance as a basic physical object and treat h and 
TZ (or 6(j) in the flat-slicing gauge) as a proper quantum field. Renormalization is then the 
natural solution to eliminate the ultraviolet divergences and keep the variance in position 
space finite and well-defined. Since the physically relevant quantity (power spectrum) is 
expressed in momentum space, the natural renormalization scheme to apply is the so-called 
adiabatic subtraction as it renormalizes the theory in momentum space. Adiabatic 
renormalization j^, |T^, removes the divergences present in the formal expression (13.11) by 
subtracting counterterms mode by mode in the integral (13.11) 



°° dk 

T 



167rGfc3 



(3.3) 



with Wk = k/a{t). The subtraction of the first term [IGTcGk^ /An'^a^Wk) cancels the typical 
flat space vacuum fluctuations, which are responsible for the quadratic divergence in the 
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integral (13. 2p . The additional term, proportional to (W^r^)*^^-* and which involves c? and a, 
is necessary to properly perform the renormalization in an expanding universe. It cancels 
the logarithmic divergence in (13.21) . 

However, one can legitimately ask if the momentum-space counterterms are uniquely 
fixed. In other words, can a different renormalization scheme, as the DeWitt-Schwinger 
subtraction prescription in momentum space lead to a different expression for the 
counterterms? Before going further in our analysis, let us briefly summarize the main steps 
of momentum-space renormalization (for an extensive exposition see {^). 



B. Momentum-space renormalization 



i. Adiabatic renormalization 



Let us consider a generic free scalar field in our spatially fiat cosmological metric 
obeying the field equation 



¥3 + 3H^ - VV + (^^ + ^R)^^ = 
where R is the four-dimensional scalar curvature 



R = 6 



+ 



(3.4) 



(3.5) 



and ^ is the curvature coupling. We include this term to be general and for illustrative 
purposes. The Fourier components of if obey the equation 



(Pk + 3-ipk + a k ipk + 
a 



6 



a2 



+ 



iPk = . 



(3.6) 



Note that the role of the dimensionless parameter ^ here is similar to the role of the e and 
?7 parameters in the wave equation (I2.29P for the perturbation. 

Let us assume that the form of the normalized modes is (2\/a^(t))~^/^e*'^^(^fc(t), where 
V = is the volume of a box of coordinate length L. The continuous limit is obtained by 
replacing V by (27r)^ and J2k / "^^^ function <^fc(t) obeys the equation 



where 



with Wk = ^y^JaP~+^rn? and 



^ 4 



a) n«^-2j^- 



~ w, 



(3.7) 



(3.8) 



(3.9) 



-1/2 if yj^(t')dt' 



The function ipk{t) is assumed to obey also the adiabatic condition i^k{t) 
in the large k regime. This condition does not uniquely determine the form of (pk{t) (different 
solutions lead to different sets of modes and, therefore, to different vacuum states) but it 
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uniquely determines an asymptotic expansion for all possible solutions. This expansion is 
characterized by 

Mt)-W~'/'e-'-f'^^^''^'''' , (3.10) 



with a recursive expansion 
where 



Wk = Wk + wf'>+wf + ... (3.11) 
-(^-.^/^^-.-^/^^-.-V) . (3.12) 



Similar expressions for w^j^^ and all the other higher-order terms can be found. Each 
order in the (adiabatic) expansion is characterized by the number of time derivatives 
of a{t) appearing in a term. While Wk is of zero order, w^^ is of second order, as one 
can trivially verify by simple counting. Note that a is of second adiabatic order. One 
should have in mind that the adiabatic expansion is an asymptotic series and does not 
converge in general. This is why the form of the modes is not uniquely defined by the adi- 
abatic condition, thus allowing to have different solutions and hence different vacuum states. 

The evaluation of the variance of the field as a sum in modes leads to a divergent 
expression. In the continuous limit it is given by 

POO 

{if") = (47rV)-^ / dkk''\^k? . (3.13) 

As is evident from the asymptotic expansion fl3.10p . one necessarily encounters ultraviolet 
divergences in the above quantity. Generically one encounters quadratic and logarithmic 
divergences. In the adiabatic renormalization, the physically relevant finite expression is 
obtained from the formal one by subtracting mode by mode each term in the adiabatic 
expansion of the integrand that contains at least one ultraviolet divergent part for arbitrary 
values of the parameters (m and ^) of the theory |j9,]. Applying this to the particular case 
of the variance of the field one gets 



oo 



{^^)ren = (47rV)-l / dkk' ' W,' - (W^'r) , (3.14) 

Jo 

where wj^^ and {W^-^Y^^ are the zeroth and second order terms, respectively, in the adiabatic 
expansion of WjT^- Generically the second adiabatic counterterm (W'^)^"^^ is given by 



^w^i — ^ 

k 



1 -1/2 d?' _i/2 1 _i 



1/2 " -1/2 , 



(3.15) 

WJ, \_Z, Ui- z, 

3ace limit (fk goes to the usual Minkowski form of the modes, {WfT^Y'^^ 
goes to zero, and {(p'^)ren gives a vanishing result, in agreement with normal ordering. In 
a dynamical universe, (Wf7^)^'^^ is generically nonzero and the renormalization produces a 
non-trivial, well-defined result for the physical variance. 

For its physical relevance in our reevaluation of the inflationary power spectra, it is 
specially interesting to analyze the massless limit of the above result. When m ^ 0, the 
subtraction counterterms take the form 

(47rV)-^A:^ + (W.'^r) = (4vrV)-^fc^ | ^ + (1 - 60;"' ^ ' 



2A;3 



a 

+ - 

a 



(3.16) 
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2. Comparison with the momentum- space representation of the De Witt-Schwinger subtraction 

terms 



Using a local momentum representation in a normal neighborhood admitting Riemann 
normal coordinates, one can express in momentum space the DeWitt-Schwinger proper time 
representation of the Green functions fl7!| (see also jl, [3]). Up to second adiabatic order, 
we have 

G'SsiP, Q) = -^Sr / ^^fce-nyrr^ + , (3.17) 



(27r)4 J ^(A;2 + m2) (k^ + m^)^' ' 

where g{P) is the determinant of the metric at the point P in the so-called Riemann normal 
coordinates y'^, ky = —k^y^ + ky, and the contour of integration in the fco plane is assumed 
to be the usual contour defining the particular two-point function. The point Q has been 
taken as the reference point for constructing the Riemann coordinates. These coordinates 
are constructed by considering the unique geodesic that joins the reference point with an 
arbitrary point P in a normal neighborhood of Q. The Riemann coordinates y^ of P are 
given by 



(3.18) 



where A is the value at P of an affine parameter of the geodesic joining Q, at A = 0, to P. 
The vector is the tangent to the geodesic at the point Q 



dX 



(3.19) 



Although in these coordinates the form of the geodesic equations is trivial, the form of the 
metric at Q is not necessarily Minkowskian. When we impose, additionally, that ds'^\Q = 
— {dy^Y + {dy^Y + {dy'^Y + {dy^Y^ we have then the so-called Riemann normal coordinates. 
In our spatially flat universe ds^ = —dt^ + a'^{t)daP' and for our problem, we take the initial 
and final points {Q = (to, x'),P = (to, x)) at the constant time hypersurface t = to. We then 
have y^ = and 

f =a(to)(x-f) . (3.20) 

Moreover g{P) = g{Q) = 1. 

The explicit form of G^^^g{P,Q) in our cosmological scenario reduces, after performing 
the ko integration, to 



G^Dsi^^ ^o; x',to) — Ii 



where R is the scalar curvature and 



h 



2(27r)3 
1 

2(27r)3 



d^k- 



d^k 



Wka 



(3.21) 

(3.22) 
(3.23) 



with wl^ 



fc2 -)- m^. Finally, performing the change of variables k A;/a(to) we get 



^Dsl^' ^oi to) 



2(27r)3a3 



d^k 



1 I il-QR ' 

Wk 2wl 



(3.24) 
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In the DeWitt-Schwinger point-splitting framework, the renormahzation of the variance 
to)) proceeds by subtracting G^^(x, to; a;', to) to the two-point function and taking 
the coincidence-point hmit x' ^ x 



(v?^(f,to))ren = limJ(v9(f,to)v?(x',to)) - G''§s{x,to]x\tQ)\ . 



(3.25) 



In the massless limit, m ^ 0, the integrand of the momentum-space representation of 
G%\{x,to;x',to) is 



^Dsl-^' ^o! ^\ to 



1 



2(27r)3a3 
1 



dkk' 



a 

2P 

a' 



+ 



+ 



sin 


k 


X — x'\ 


k 


X 


— x' 





Taking x —>■ x', the integrand of the above expression coincides exactly with the momentum- 
space counterterms obtained previously via the adiabatic renormahzation. Therefore, the 
question ending section IIII Al is answered. The adiabatic subtraction coincides with the 
DeWitt-Schwinger subtraction prescription in momentum space in the massless limit, irre- 
spective of the value of the dimensionless parameter ^. 



C. Tensorial spectrum 

Let us apply the above scheme to the tensorial fluctuations of the metric. To determine 
the counterterm (VF^^)^^-* one should rewrite the wave equation fl2.3p in the form (13.71) . This 
can be easily obtained by performing the change = a^/'^h^. One obtains = k/a (i.e., 
m = 0) and a second order adiabatic term a of the form (13.91) with ^ = 

3 / d\ ^ 3a 



A straightforward calculation gives 



a2 



iW,TT^ = ^, + ^,. (3.28) 



To perform the explicit computation it is useful to take into account the slow-roll relations 

da (1 + 2e) 



dr {Ht'^) 
d^a (2 + 5e) 



dr2 (Ht^) 

These derivatives are related with the usual dotted derivatives as follows: 

_i da 
a = a — 
dr 



(3.29) 
(3.30) 

(3.31) 



-^''hz) +^'"^. (3-32) 



Kdr I dr 
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The final expression for zerotli and second order adiabatic counterterms is tlien^ 

(2 + 3e)" 



w 



-1 



(3.33) 



Tlie subtraction of these counterterms produces the expression (assuming that E{k) = 1 
and F{k) = for every k)^ 



dk IQTiGk^i-TTi) 



k 



Hi'\-kT)\' 



-K^—kr) 



1 + 



(2 + 3^ 
2kH'^ 



Therefore, the renormalized expression for A\{k^t) is 



Anna? 



m\-krf 



1 + 



(2 + 3e) 
2Pr2 



7r(— A;r) 

Note that the asymptotic behavior for large k of the Hankel function is 



|if«(-A:r)p 



7r(— fcr) 



1 + 



(4z/2 



1^ 



8(-fcr)2 _ 



+ ... . 



(3.34) 



(3.35) 



(3.36) 



The first two terms for z/ = 3/2 + e are exactly the adiabatic counterterms subtracted above. 

It is useful to remark that, in adiabatic renormalization (or, equivalently, in the DeWitt- 
Schwinger subtraction algorithm), the number of terms that should be subtracted is deter- 
mined by the degree of divergence of the expectation value under consideration. In our case, 
this degree of divergence is of second-order. Therefore, one must stop necessarily at second 
order if the purpose is to make finite the expectation value {h'^)- The relevant subtractions 
are carried out mode-by-mode and result in a sum-over-modes that has no ultra-violet di- 
vergence. In addition, as we have already stressed, it is important to point out that it is 
necessary to subtract the relevant adiabatic counterterms for all modes, including those with 
large wavelength, even if the expansion of the universe is not slow. Note that the reason for 
the name "adiabatic" regularization is a result of the fact that the adiabatic (or arbitrarily 
slow) limit is used to identify the necessary subtractions. Therefore, it is not correct to 
argue that the adiabatic renormalization terms are invalid when the expansion is rapid, i.e., 
in the non-adiabatic regime. 



D. Scalar spectrum 

We will work out the (renormalized) spectrum of scalar perturbations using the gauge- 
invariant formalism. To convert equation 02.371) into one of the form 03.71) we have to 



2 



In a similar way, the DeWitt-Schwinger subtraction terms in momentum-space lead to the same result. 



Note in passing that neither of the proposed (adiabatic) counterterms defined in [28| for the tensorial and 
scalar power spectra seem to agree with those obtained here by adiabatic regularization. 
^ Note that the second adiabatic counterterm introduces an infrared divergence in the continuum limit. For 
fc — > the integrand in {h'^)ren is proportional to dk/k. This logarithmic divergence disappears in the 
finite box formulation and has no impact on the observable power spectrum, in contrast to the ultraviolet 
divergences. Furthermore, the calculation of the renormalized stress-energy tensor (T^^), and hence the 
evolution of the universe, is insensitive to the infrared cutoff L when L — > oo. 
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perform the change TZk = a^^'^zlZk- We then obtain that Wk 
order adiabatic function a is given by 



fc/a, and that the second 



4 \a 



(3.37) 



Note that this a is different from that of fl3.9p . The reason for this can be easily seen by 
looking at the wave equation in the flat-slicing gauge (12.291) . That equation does not have a 
constant massive term, and the second adiabatic order term is if^(3?7 — 6e), which is just the 
last term in (13.371) . This implies that the scalar perturbation (the inflaton in the flat-slicing 
gauge or IZ in the gauge- invariant approach) should be regarded as a massless field. ^ The 
corresponding second order counterterm is thus of the form 



+ 



2awl 



Zwl a 



(3.38) 



Therefore, the final expression for the counterterms is 



1+ 



(2 + 3(3e-77)) 
2(-fcr)2 



(3.39) 



Proceeding in a parallel way as for the tensorial case, we get a renormalized value for 



k 



(3.40) 



Therefore, 

ren ~ 



°° dk , o -Trr 
k 4(27r)^z2 



Tl{ — kT) 



1 + 



(2 + 3(36-7/)) 
2(-A;r)2 



(3.41) 



where the index yU is given byyU = 3/2 + 3e — 77. The renormalized expression for A^(/c, t) is 



A7Tk^{—'n-T) 
4(27r)3;z2 



Hi'\-rk)\'- 



Tl{~kT) 



1 + 



(2 + 3(3e-r/)) 
2{-kTV 



(3.42) 



Note finally that, as expected, the same result is obtained working in the fiat-slicing gauge. 
One then obtains 



fc3(— ttt) 
Anna? 



T^i—kr) 



(2 + 3(36-7/)) 

2{-kTf 



(3.43) 



where /i = 3/2 + 3e — 77. Taking into account that A^(/c,t) = {0? / z^)^^^{k,t) one recovers 
the expression fl3.42p for A^(fc, t). 



^ This was first taken into account in the calculation of the renormalized power spectra in 0, 0| in terms of 
the gauge-invariant quantity TZ. This quantity was also quantized in [30|, but there the V" term in (|2.27p 
was treated as zeroth adiabatic order, unlike the present treatment, in the evaluation of the counterterms. 



18 



IV. TESTABLE CONSEQUENCES 

In this section we work out new expressions for the observable magnitudes that follow 
from the renormalized power spectra of the previous section. 



A. Tensor-to-scalar ratio 

Let us now consider the tensor-to-scalar ratio 



4^ 



4 



WirGk^-Tir) 4(27r)V 




"l 1 (2+3e) " 

^ 2{-fcT)2 


) 




47r22a2 47rA;3(-7rr) ^ 




\. (2+3(3e-r,))" 

^ 2(-fcT)2 


) 



Since at the Hubble exit time tk we have z'^{tk) = 2Mpe(tfc)a^(tfc), it follows that 



4- 



t—tk 



We also find that 



\Hl'\-rk)\ 



(2 + 3e) 



16e(tfc) 



7r(— fcr) 



1 + 



t=tk 



IT 



'-ae{tk 



(4.1) 
(4.2) 

(4.3) 



where a is a numerical constant of order unity, a ~ 0.904. A similar estimation is obtained 
for the factor coming from the scalar spectrum (up to the substitution — > /x, i.e., e — >■ 3e— ry) 



H;!\-rk)\'- 



7r(— A;r) 



1 + 



(2 + 3(36-7/)) 
2(-fcr)2 



t=tk 



-aiSeitk) - vitk)) ■ (4.4) 

TT 



Therefore, the tensor-to-scalar ratio r, evaluated at the Hubble exit time, is 



16e(4) 



e(4) 



3e(tfc) - ?7(tfc) 



(4.5) 



An important comment is now in order. To obtain this result we have evaluated the 
power spectra at the Hubble exit time. Nevertheless, as we will show in the next paragraph, 
this estimate does not depend critically on the precise time at which the counterterms are 
evaluated. During slow- roll inflation the counterterms decay as Irp*^ (tensorial), jrp^^'^"'') 
(scalar), in the late-time limit as |r| — > 0. This means that they decay very slowly and 
are in fact constant in exact de Sitter inflation. After reheating, the counterterms decay 
more rapidly, as is obvious from fl3.16p . One would then recover the standard prediction 
r = 16e(tfc). However, since the modes acquire classical properties soon after exiting the 
Hubble sphere, the relevant time to evaluate these magnitudes falls in the interval between 
tk and a few e-foldings after it. 

Let us analyze in detail the time dependence of A^(A;,t), A'^{k,t), and r in terms of the 
number n of e-folds after crossing the Hubble radius. Since —kT{tk) = 1 + e, we can write 
—kr as 

—— = (1 + e)(^^) = (1 + e)e . (4.6) 
T{tk) a 



-kT = (l + e] 
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The expression for Al{k,t) when —kr ^ 1 leads to 

2 



[1 - ^-kr^'^^r^] > (4-7) 



and using (14.61) . we find 



Ml\ 2ti J ^ ^ ' 2(-/cr) 



^S(*-")-Mi^Vll-^e-T (4.8) 



Assuming that we are just a few e-folds after the Hubble exit but before the end of inflation, 
i.e., n > 1 but ne -C 1, so that e"^*^" 1 — 2e?2, we obtain^ 



K{Kn) ^ j e(4)(2n-3/2) . (4.9) 

A similar estimation can be obtained for the scalar power spectrum 

Note that the parameter n enters in the power spectra parameterizing the (unknown) time 
at which the modes exhibit classical behavior. However, since both tensorial and scalar 
spectra have the same dependence on n, the tensor-to-scalar ratio is not sensitive to the 
unknown parameter n, in the same way as it is insensitive to the scale of inflation H{tk), 
and it is essentially given by (14.51) . The same conclusion can be drawn if one estimates the 
time derivative of r. One then obtains that 

|-.0(.V (4.11) 

This means that the renormalized value for r is changing very slowly, in agreement with 
the previous evaluation. In other words, the value of r evaluated at the Hubble radius 
crossing time remains nearly constant during this period of inflation. Obviously it is not 
strictly constant, as in the computation without renormalization, but its change is slow.^ 
We regard this result as a signal of the robustness of renormalization in determining the 
spectra of inflationary perturbations. 



^ We note that if one evaluates the power spectra at the end of the slow-roll era (where ne ~ 1) the 
contribution of the counterterms is still significant. However, we find it more natural to evaluate the 
spectra soon after Hubble exit, when the modes have already acquired classical properties. 

^ To be precise, d{\nr)/dn — 2e — rj. Note that for the exponential potential model V{(j)) oc 
exp[— ((/)/Mp)a/2/p] we have r/ = 2e. Then the new predictions coincide, by accident, with the stan- 
dard ones. So in this case f = 0. 
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B. Spectral indices 

The above calculations have another important consequence: the spectral indices remain 
unchanged when they are evaluated at tk or a few Hubble times after it. By definition, we 
have 

^ ^ d\nPt ^ ^ d\nH{tk) _^ (ilnefa) 
dink dink dink ' 

which, according to (14. 9p . turns out to be independent of n. The result is 

nt = 2{e-T]), (4.13) 

where e and r] are evaluated at t^. For the scalar index we have 

n -1= "^^^^^ = o ^lngfa) _ dlnejtk) ^ ciln(3e(tfc) - r/fa)) 
In /c d In /c In /c d In A; 



Explicit evaluation gives 



(12e2 - Sen + n , , 

n,-l = -6e + 2T] + ^ '- — ^ 4.15 

3e — r] 



where ^ is another slow roll parameter lOl]: ^ = Mp{V'V"' /V"^). This parameter can be 
reexpressed in terms of e,?7 and the running of the tensor ial index = drit/dlnk as follows 

n[ = 8e{e-r]) + 2^ . (4.16) 



C. Consistency condition 

The formulas (14.5^ 14. 13114. 15]) provide an algebraic relation between the tensor-to-scalar 
ratio and the spectral indices, r = r{nt,ns,n[), that takes the following form 

r = 4(1 - - fit) + ^at'L ~ ~ V^^^ + ~ ~ " ^^"^'^^ 

Positivity of the argument of the square root in this expression imposes the following con- 
straint 

{l-Usf -n^t+2n[>0 . (4.18) 

The consistency relation fl4.17p . contrasts with that obtained without invoking renormal- 
ization, namely, r = —8nt. Since r must be positive, according to the standard prediction 
one would expect rit to be negative. This restriction, however, does not follow either from 
f l4.17p nor from fl4.18p . Even more, from (14.180 it is easy to see that our prediction allows 
for Tit > 0, in sharp contrast with the standard prediction, which necessarily requires 
Ut = —2e < 0. We also remark that, according to the standard derivation, the running of rit 
is fully determined by the values of and Ut, since then one finds n'^ = — nj(l — -t-n^). On 
the contrary, the manipulations that led to (I4.17P indicate that n'^ is now an independent 
quantity that needs to be measured in order to check the new consistency relation (14.171) . 
This aspect could make more challenging the experimental verification of the consistency 
condition of (single-field) slow roll inflation. 
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We stress again that the new consistency relation allows for a null tensorial tilt Ut = 
while being compatible with a non-zero ratio r ^ 4(1 — n^). Since the observations from 
the 5-year WMAP [3i](with BAO+SN) strongly suggest that (1 - n,) ^ 0.030±0.015 (with 
r < 0.22)^ it follows that, for the case of an exact scale invariant tensorial power spectrum, 
rit = 0, our consistency relation leads to r ^ 0.12 ± 0.06. These possibilities may soon 
come within the measurement range of forthcoming CMB polarization experiments such as 
the PLANCK satellite, the QUIJOTE CMB experiment, SPIDER, Polar BEAR, EBEX, 



BICEP, SPUD, and the future CMBPol mission [31 



D. Comparison between the standard and the new predictions 

We can compare the differences between the standard and the new predictions in terms 
of the representative set of chaotic potential models V{(f)) = XM^p'^^cff. The standard 
prediction is 



r-^ l_n=i^±^ n, = -^ n' = ^ (4 19) 



where N = lnaend/ct{tk) is the number of e-folds of inflation between the horizon crossing 
time tk of cosmological wavelengths and the end of inflation. If we invoke renormalization 
we get, instead, 

- '"^ 1-".-^ . ".-^ . (4.20) 



{p + 2)N ' 2N ' 2N ' ' 2N^ 

Note that now the quadratic potential is characterized by having an exact scale-invariant 
behavior for the tensorial spectrum. 



E. Comparison with WMAP data 

Although a definite test of the new predictions, like the above proposed consistency 
condition, requires more accurate data on tensor fluctuations (so far we have only upper 
limits on r) we can compare the new predictions with the standard ones on the basis of 
the five year WMAP data. We can contrast the predictions (14.201) for r and Ug, for the 
representative values p = 2 and p = 4, with the WMAP 5-year data (see Fig. 5 of Ref. j^). 
We find that both models are compatible with the experimental data for the reasonable 
range of N between 50 and 60. This is in sharp contrast with the prediction of the standard 
approach, where the monomial potential with p = 4 is excluded convincingly. 



V. CONCLUSIONS AND FINAL COMMENTS 



Inflationary cosmology predicts that, due to quantum effects, small density perturba- 
tions are generated in the very early universe with a nearly "scale-free" spectrum. The 
detection and analysis of anisotropies in the cosmic microwave background has confirmed 



^ Marginalized over all other parameters of a flat ACDM model. 
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this prediction. Moreover, inflation also predicts the creation of primordial gravitational 
waves, which still remain undetectable. Forthcoming high-precision measurements of the 
cosmic microwave background [31j may measure effects of relic gravitational waves, and this 
will be crucial to test the inflationary paradigm and strongly constrain inflationary models. 
Therefore, it is of crucial importance to scrutinize, from all points of view, the quantitative 
predictions of inflation. 

In this work we have pointed out that, if quantum field renormalization is taken 
into account, the predictions of (single- field) slow-roll inflation for both the scalar and 
tensorial power spectra change signiflcantly. In our physical context, renormalization is 
naturally implemented in a mode-by-mode subtraction scheme and uniquely deflnes the 
momentum-space counterterms. These counterterms are evaluated in the period when the 
perturbations acquire classical properties and this leads to testable predictions that differ 
signiflcantly from the standard ones. Because of the question of the underlying quantum 
nature of the gravitational fleld, there are several ways one could think about the process 
by which the dispersion spectrum of the inflaton fleld influences the spectrum of the scalar 
perturbations of the metric. One way is to think of the gravitational fleld as making a 
measurement of each mode of the inflaton fluctuation fleld within a few Hubble times of 
its exit from the Hubble sphere. If this measurement of the inflaton dispersion spectrum 
is similar to standard measurements in quantum mechanics, then it should measure the 
renormalized value of the inflaton spectrum at the time when the measurement is carried 
out, which we take to be the time when the perturbations acquire signiflcant classical 
properties. The time t at which this occurs affects both the scalar and tensorial power 
spectra, but for a large range of values of t, the effect of renormalization remains signiflcant 
for wavelengths that today are at observable scales. We have shown that if the power 
spectra are evaluated n e-folds after the Hubble radius exit time, where n -C 1/e, then the 
observable parameters r and Ug, rit, n[ are insensitive to the value of n. If the power spectra 
were evaluated at a time t well after the end of inflation, one would recover the standard 
predictions. With the present understanding of the nonlinear aspects of quantum gravity, 
it is difficult to reach a definitive answer regarding the value of ra, so the fact that we find 
observable differences offers a deep way to experimentally probe this question. 
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